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ABSTRACT 


Design, fabrication and cold flow testing of a modeled 
jet engine test facility was conducted in an effort to 
provide an inexpensive vehicle to study geometric variations 
in diffuser geometry which could improve system efficiency. 
The design is based on Mach number similitude and consists 
of two configurations currently in use at the Naval Air 
Propulsion Center, Trenton, New Jersey. A constant area 
diffuser and a variable area diffuser with translating 
centerbody were modeled. Baseline mapping of the operating 
characteristics for each diffuser with representative scaled 
engines was conducted to provide a reference against which 
alternative geometries would be evaluated. The constant 
area plus two variants were tested. A five-sixths and two- 
thirds reduction were studied to investigate the potential 
femeinecreasing efficiency for a specific engine diffuser 
Combination at NAPC. Secondary flow provisions were incor- 
porated into the design to allow variation of this parameter. 
The modeling results were consistent with theory and the 
test apparatus produced repeatable results. A two dimensional 
double ramp (wedge) capable of being translated in a rectangu- 


lar duct was suggested as an alternative diffuser geometry. 
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I. INTRODUCTION 


ie by LO elticiently exereise control over the energies 
entrained within a supersonic airstream has been the quest of 
aerodynamicists for several decades. The designers of wind 
tunnels, jet inlets, gas dynamic lasers and jet engine test 
facilities have each addressed the gas dynamics of this topic. 
Hach design has had to incorporate a method to decelerate the 
flow, generally, through a mechanical device such as a diffuser. 
The complexities of treating the recompression of a real fluid 
in the presence of a boundary layer have defied analytic 
modeling of a supersonic diffuser to any great extent. The 
design approaches taken have been empirically based, which has 
led to a wide variety of diffusers tailored to meet the unique 
Operating environment at a particular facility. This study 
is sponsored by one such facility challenged with one of the 
consistent fascinations of modern engineering: how to extend 
the limits of one's design in the presence of new technology 


Or shifting economic variables. 


A. PROBLEM FORMULATION 

Ground testing of jet engines has long been an integral 
Part of the design and maintenance practice in both the military 
and commercial avaition industries. Organizations, chartered 
with the testing of these engines, strive to generate a test 


envelope which closely approximates the operating envelope 


To 





which the engine will encounter in service. Advances in engine 
technology have imposed added demands upon the test engineer to 
extend the test envelope accordingly. This challenge has 

proven a classic cost effectiveness exercise, wherein, as 

higher altitude testing at increased power is pursued from 

one end of the spectrum, the attendant cost of exhausting the 
effluent in an innocuous manner to the environment spirals. 

The economic challenge continues to compound over the life cycle 
of the facility as energy costs associated with demands on 

the exhausters escalate. 

Test cell philosophy has focused foremost on achieving a 
sufficientiy flexible design which will accommodate a wide 
range of engines. Large exhaust mechanisms, capable of 
handling a wide range of exhaust states, were adequate when 
the motive energy cost was only a small fractional cost of 
the total price of testing. Strategies to enhance pressure 
recovery prior to the exhauster were developed but optimiza- 
@eeomeert the design in this regard was not a bonafide concern. 
The present testing scenario reveals that the associated 
costs in exhausting the effluent rivals any of the other cost 
variables and percentage improvements in efficient 
pressure recovery through retrofit of the original design 
merit consideration. 

A typical test cell design is as depicted in Figure i. 

The engine to be tested is mounted on a test bed and located 
in the test cell such that the exhaust will be vented into 
an augmenting tube which acts as an ejector-diffuser assembly. 


la 





The kinetic energy of the exhaust stream is converted by the 
diffuser into a pressure for presentation to the exhauster. 
Each cell is nominally equipped for secondary flow in which 
secondary air is entrained with exhaust jet gas to provide 
engine cooling and dilute the combustion products. Allowance 
is made for relative positioning of the test bed and diffuser 
to reconcile potential problems with pressure gradients under 
conditions of secondary flow which may influence the operating 


envelope. 


B. LITERATURE SURVEY 

Several searches were conducted to survey the available 
literature for supersonic ejector-diffuser studies and 
Theoretical discussions germaine to this investigation. An 
online computer search of several national data bases was 
conducted using the keyword, keyphrase approach. Results of 
the search revealed over 10,000 documents generally associated 
with the broad topic area,of which,a highly focused search 
indicated over 300 documents with relevant material. A 
hardcopy of the latter with a brief synopsis of each report, 
was procured for further review. The survey was restricted 
to English or English translations but evidence of many 
feGeiegn papers on the subject was apparent. In no respect 
is the review considered all-inclusive. 

A synopsis of the most recognized works gives a flavor 
ferethe approach adopted. In 1949, pioneer work, which 


appears as a baseline in most studies related to supersonic 


ae 





diffusers is attributed to Neumann and Lustwerk {Ref. 1}. 
This study included a one-dimensional theoretical analysis, 
and an experimental modeling with flow visualization by 
Schlerin photography of a constant area diffuser. A "tran- 
sverse shock" was observed and categorized as the operative 
Meenanism controlling diffusion. An optimum diffuser with 
an L/D of 10 was identified. In 1958, Lukasiewicz {Ref. 2}, 
studied data from several existing wind tunnel diffusers, 
concluding that fixed geometry diffusers can approach the 
pressure recoveries established from normal shock theory. 
Pressure recovery, far in excess of that obtainable with a 
constant area diffuser, was established for systems which 
employed variable area diffusers. In 1954, Hastings (Ref. 3} 
established the beneficial effect in diffuser performance 
of auxiliary ejection to partially evacuate test cells. 
Numerous additional studies with specific design goals have 
been conducted to optimize test facility operation. The most 
extensive noted were those conducted by Panesci and German 
{Ref. 4}, for Arnold Engineering Development Center in the 
60's in which variable geometry diffusers with a centerbody 
Memes cWployed. Here again, pressure recovery far in excess 
of that achievable with constant wea diffusers was observed. 
Generalized studies to characterize pertinent parameters 
governing the flow phenomenon in rectangular diffusers were 
conducted by Merkli {Ref. 5} in 1976 and Waltrip and Billig 
{Ref. 6} in 1973. Merkli focused upon Mach number, diffuser 
length, boundary layer and Reynolds number as controlling 
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parameters. Reynolds effects were discounted as minimal with 
Mach number and diffuser length the significant parameters. 
Waltrip and Billig corroborated previous works establishing 

8 - 12 tube diameters as the required recovery zone. They 
also focused on an oblique shock system as the governing 
mechanism. 

Ginoux {Ref. 7} compiled an excellent summary of a short 
course in Supersonic Ejectors conducted at the von Karman 
Institute. The short course was an attempt to focus on the 
most advanced initiatives and progress in theoretical modeling 


and design of high efficiency ejectors. 


GC. RESEARCH OBJECTIVES 

The Naval Air Propulsion Center at Trenton, New Jersey, 
as a major jet engine test facility, has experienced the 
technological advances in engine design which have approached 
the design limits of their ejector-diffuser assemblies 
exacting a heavy burden on power consuming exhauster machinery 
to maintain simulated altitudes. As an adjunct to a much 
Hemeer study, a cold flow modeling of their existing plant 
was sponsored by the center. The principal goal, assuming 
Satisfactory modeling of the test facility, vas to test 
alternative diffuser geometries in anticipation of enhancing 
overall efficiency. The modeling process was such that, 
Mach number similitude could be maintained, any efficiency 
increases over the baseline would be due largely to Poona tne 


effects. 
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As discussed in the general treatment on diffuser theory, 
the two principle types of diffusers, a fixed area and a 
variable area, were modeled. In both diffusers recompression 
of the supersonic flow is accomplished by a complex-shock 
mechanism under the influence of a boundary layer, with post 
shock subsonic diffusion following recognizable theory. The 
experimental technique devised was to establish the diffuser 
characteristic on a non-dimensional basis as a baseline 
against which ‘new! geometries may be judged. Operating 
envelopes for each diffuser design would be duplicated as far 
as practicable with the same engines. Whereas the phenomenon 
by which recompression occurs would not be directly studied, 

a pressure histogram along the diffuser was recorded in order 
to postulate the character of the operative mechanism. It 

was anticipated that attempting to control the shock mechanism 
would likely provide the largest gains in efficiency as opposed 
to manipulating the subsonic ieee cen process. 

The scope of the investigation would be guided by studying 
only those configuration changes which could readily be 
retrofitted into the existing space limitations of the parent 
facility. Conceptual designs would be unbounded by any 
environmental or stress-related constraints, allowing a 
sponsor's cost benefit analysis to sort out those aspects 
of new design proposals. 

Despite successful construction of a highly flexible 
model, a major portion of the stated objectives could not 


be accomplished within the timeframe alloted to this phase of 
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the study. As baseline testing proceeded into the variable 
geometry diffuser, Figure 2, an unanticipated heating and 
vibration phenomenon was observed. The extent and nature of 

the phenomenon was not readily ascertainable but was in evidence 
only with the use of the centerbody configuration. The problem 
was of such proportion as to potentially taint the conclusive- 
ness of future work involving devices imbedded in the jet 
stream. A separate detailed study of the phenomenon was 

ordered and a new set of objectives was established in concert 
Mies the sponsor. 

in an effort to optimize test cell geometry for one of the 
more heavily tested engines (F404), a series of liners which 
would reduce the cross section for diffusion were designed 
for insertion into the full scale straight tube diffuser. 
Engines were tested in anticipation of achieving better diffuser 
efficiency by seeking to optimize the ratio of nozzle area to 
Seeniicer cross section for the highest pressure recovery. 

The details of the model design and testing in the context 
of this narrower objective are contained in the thesis proper. 
The conceptual work related to the original objective with 
a proposal for an alternative method of diffusion are discussed 


in Appendicies A and C. 
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Pte eb eORY AND ANALYSIS 


Pressure recovery in a supersonic jet engine test facility 
is accomplished by a mechanical device called a diffuser. Two 
types of diffusers are recognized, the fixed or constant area 
diffuser and the variable geometry diffuser. The fixed geometry 
normally is associated with fairly constant input parameters 
such as mass flow rate, stagnation temperature and pressure. 
The variable area geometry is utilized where fluctuations in 
fluid characteristics or engine geometry (such as variation in 
exhaust area accompanying jet engine testing from the non- 
afterburning to afterburning mode) are an integral part of the 
testing. Hach type of diffuser may serve an ancillary role to 
eject secondary air used in cooling the engine assembly and 
tees cell. 

In each diffuser the operative mechanism which accomplishes 
i@emtirat order pressure recovery from supersonic to subsonic 
conditions is a shock system. Subsequent pressure recovery must 
morimow the guidelines for subsonic diffusion. Projecting an 
improvement in efficiency accompanying any alternative geometry 
would require a projection of the probably shock patterns and 
the interaction of that shock system with a postulated boundary 
iijer., ihis interaction, in simple geometries, has not been 
conclusively researched; hence, this type of approach in the 
presence of complex geometries is not warranted. Analytic 


models to guide the design of a new geometry for jet engine 


lee 





testing abound in the literature but generally assume the most 
convenient of assumptions. The model is generally one dimen- 
Sional steady state using a simplified control volume and 
serves to bound the expectations only. Academic interests 
aside, a purely empirical approach is warranted. The approach 
adopted herein calls for establishment of baseline models of 
proper similitude with the existing facility from which 
Characteristic curves can be drawn and against which alter- 
native designs may be mapped and contrasted. 

Acceptance of any observed change in system efficiency 
Merics consideration only if dynamic similarity of the flow 
field has been verified between the baseline model and the 
emens facility. With supersonic compressible flow, Shapiro 
(er. S!, is replete with support documentation illustrating 
the role of Mach number as the significant parameter in 
@feracterization of the flow. Merkli {Ref. 5}, in a series 
of experiments with rectangular constant area supersonic 
diffusers, concluded that Reynolds number has little effect 
on the pressure recovery. Mach number, as the ratio of kinetic 
energy to internal energy, was thus chosen as the best para- 
meter upon which to base model development. Geometric 
compatibility was governed by the constraints of the engines 
to be tested and the limits imposed by the available air supply 
at the model test facility. The influence of temperature 
between cold ambient testing and prototype testing with hot 
exhaust gases would be addressed in the discussion of results 


as how it might impact the operative pressure recovery 


18 





mechanism. Appendix A provides a more detailed study of 


modeling/scaling considerations peculiar to this study. 


A. FLOW CHARACTERIZATION 

Flow at the exit plane of the nozzle achieves supersonic 
proportions whose Mach number may be approximated by analyzing 
a Prandtl-Meyer corner flow from the nozzle exit to diffuser 
entrance. The increase in area from the exit plane to the 
diffuser allows the jet to expand supersonically as it fills 
the available volume. A Prandtl-Meyer expansion may also 
be utilized to estimate the pre-shock Mach number. Shocking, 
due to perturbation of the jet stream with the boundary layer 
as reported in several other workd {Ref 8} will assume an 
oblique character. The oblique shock system will, upon 
eeenainment Of subsonic conditions, blend into a turbulent, well- 
mixed stream which would diffuse in accordance with subsonic 
theory. The oblique shock system would be expected to 
migrate along the diffusers length for a given geometry of 
diffuser, in some proportion to the driving pressure. The 
oblique shock system, as discussed by Shapiro, {Ref. 8}, 
will either be strong or weak as governed by the stability of 
the flow, the nature of the boundary layer interaction and 
a multiplicity of lesser related factors. Pressure variations 
caused by area change conceivably promote an alternating 
compression and expansion character to the flow wherein the 
jet may tend to pulse. Restricting the flow to a constant 
area would tend to damp out this type of behavior. Figure 3 


illustrates the anticipated character of the flow. 
ilies 





B. OPERATING CHARACTERISTICS 

The operation of the variable area ejector-diffuser 
provides insight into the complexities involved when designing 
Or redesigning a new pressure recovery mechanism. Utilizing 
the simplified arrangement of Figure 4 to guide the discussion, 
the operation of this device may be described. As the total 
pressure is increased, flow in the nozzle accelerates until 
sonic conditions (M=1) are attained at the throat. Increasing 
total pressure or holding total pressure at this level and 
reducing the back pressure will cause a normal shock to stand 
in the supersonic region of the nozzle. A further lowering 
will cause the shock to pass into the test cell and into the 
diffuser. With a second throat, once the shock has been 
swallowed, the diffuser is considered started after which 
exhaust pressure may be raised shifting the shock to gones 
Where stagnation pressure loss is less. A minimum loss will 
Mreeur if the shock is located at the second throat. This 
may be accomplished by adjusting the axial position of the 
eenterbody. The minimum flow area of the diffuser, Ad, must 
be greater than A* or the cell would become choked and altitude 
Simulation could not proceed. The band of pressures, where 
cell pressure is independent of exhaust pressure, establishes 
the operating range of the diffuser. Conservatively, the 
shock is maintained upstream of the throat to preclude reverting 
moma higher cell pressure due to fluctuations in the flow 


field. 
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G. ENERGY CONSIDERATIONS 

The presence of a shock wave arising from the supersonic 
starting process represents an increase in entropy at the 
expense of stagnation pressure. The entropy rise (pressure 
loss) is greatest across a normal shock as opposed to that 
across several oblique shocks. A simple illustration using 
Figure 5 makes the point. For a flow of Mach 2.0 at the 
diffuser entrance, a one-dimensional normal shock gives a 
Stg@enacion pressure ratio across the shock of .721 with a 
post shock Mach number of .475. Using a device to diffuse 
the flow in oblique steps, then allowing for a normal shock, 
should increase the stagnation pressure ratio compared to 
the normal shock alone. 

Choosing turning angles of 6 degrees for each of two 
successive redirections of the flow followed by a gradual 
turn prompting a normal shock yields an overall stagnation 
Mmeecsuire ratio of .951. The pair of oblique shocks increases 
the stagnation pressure rise by a factor of 1.32. in the 
miomib, an infinite number of small oblique shocks will tend 


towards an ideal recompression. 
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III. EXPERIMENTAL APPARATUS 


Each of the scale model altitude test facilities 
constructed consisted of a common test cell, and an exhaust 
plenum with a variable diffuser assembly as illustrated in 
Figures 6 and 6a. Primary and secondary air were provided 
by a common source, an Allis-Chalmers twelve stage axial 
compressor. Exhaust plenum pressure was controlled by an 
air ejector driven by the common air supply from the axial 


compressor. 


A, TEST CELL/ENGINE ASSEMBLY 

The test cell, Figures 7, 7a, and 8, housed engines and 
provided a plenum for secondary air flow. The cell was 
fabricated from aluminum and of cylindrical design measuring 
[eeinehes in length and 12 inches in diameter (1.D.). The 
upstream flange assembly (1) provides a mating surface for 
igemprimary air piping, structural support for a cantile- 
vered engine housing (2) and an air seal assembly. Dry 
silicon rubber seals guarding against air intrusion are 
prescribed owing to the vacuum created for altitude simula- 
fom) A 3 inch diameter penetration (3) at the base is 
provided for secondary flow connections. The downstream 
flange (4) accommodates diffuser assembly attachment and 
incorporates a similar air sealing arrangement. Ports for 
direct sampling of cell pressure and remote connectors for 
engine pressures were provided. 
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The engine assembly, also of aluminum, consists of 3 
ici i. D.) entrance piping (5) which in addition to its 
flow straightening function served as the support for the 
engine mounting assembly (6). The mounting assembly served 
to transition the flow from the entrance piping to the 2 
inch (I.D.) conformal entrance duct. The mounting assembly 
introduced one element of versatility via a variable spacer 
Ping (7). The spacer ring allows for 2 inches of horizontal 
realignment of engines should variation in standoff distance 
to the diffuser be required. The engine mounting surface 
(8) was machined to provide a retaining collar and indented 


for set screw assembly of engines. 


B. #XHAUST PLENUM 

Interfacing beweeen the exhaust air ejector assembly and 
the diffuser assemblies was an exhaust plenum 3 foot by 3 
foot in cross section by 4 feet long. The plenum houses a 
remotely operated traversing mechanism used to dirve the 
multiple angle centerbody assembly which is peculiar to the 
variable diffuser geometry. A maintenance access/inspection 
Merc is provided to assist in alignment. A six inch access 
connects to the air ejector piping to provide closure with 


the atmosphere and a means of back pressure control. 


C. DIFFUSER ASSEMBLIES 
Two scaled diffuser assemblies, Figures 2 and 9, were 
developed to establish the baseline against which alternative 


designs can be compared. 
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ieeeotraight Tube Ditfuser 


Memnodalmceletcsus Of a |5.2> inch long 2.71, inch 
(I.D.) cylindrical ejector-diffuser. Pressure taps were 
installed to record the pressure recovery process and are 
MmeiGetrated in Figure 10. Taps were placed at one (1) inch 
intervals along the length of the diffuser. Sealing was 
achieved by rubber seals in the end flanges. The length to 
diameter ratio was 5.62. 

Two variations of this geometry, Figures 9a and 9b, 
were developed to investigate extending the operating envelope 
of the test cell to enhance efficiency and economy of operation. 
Meme@epicted in Figures 9a, 9b, and 11 inserts were added to 
achieve a 5/6 and 2/3 reduction in diameter. Two end inserts 
(9) were included to allow investigation of sudden expansion 
versus gradual diffusion in the end section. 

2. Variable Area Diffuser 

Variable area diffusion was developed by traversing 
a multiple angle conical centerbody (Figures Z and 12) within 
feeeeonech long cylindrical to conical diffuser. The overall 
length to inlet diameter ratio was 6.92. 

The centerbody was 16.5 inches long having a leading 
cone of half angle 19.8 degrees and three trailing truncated 
cones of 10.8, 8.9, and 2.6 degrees, respectively, with a 
cylindrical afterbody. Centering was provided by a reinforced 
spider (10) which provided bearing support for 3/4 inch steel 


drive shaft. The shaft was coupled to an electrically operated 
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drive mechanism, Figure 13, which was remotely activated, 
allowing travel of 6 inches with positive mechanical and 
electrical limits. Positioning circuitry generated a plus/ 
minus 5 volt output which is remotely retrieved at the 
principal operating station. 

The cylindrical to conical diffuser (8 degrees half 
angle) was equipped with static pressure taps longitudinally 
located along the wall, as shown in Figure 14. 

The integrated centerbody and diffuser permitted 
eeemrartation in the flow area presented to the jet, including 
Meemecret ion of a variable second throat. Variation in flow 
area with axial position of the centerbody is shown in 


Figure 15. 


D, &NGINES 

Two sets of engines, Figure 16, were developed to model 
the F404 and the TF30 engines tested at NAPCG. The engines 
were scaled to simulate the IRP and max A/B mode of testing. 
IRP represents Intermediate Rated Power which represents the 
highest power level without afterburner. This term is used 
Synonymously with non-afterburning throughout the thesis. 


A/B refers to the maximum afterburning mode. 


eee elk SUPPLY 
Compressed air from the Turbopropulsion Laboratory's 
Allis-Chalmers, twelve stage axial compressor, Figure 17, 


was utilized in all model testing. Maximum discharge pressure 
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of this machine was approximately 3.0 atmospheres at 15.0 
lbm/sec mass flow rate. 

Primary and secondary air, as previously shown in. 
Figure 6a, were supplied to the engine model and test cell, 
meeectiyely, through three inch I.D. piping. A six inch 
I.D. suction line was attached to the exhaust plenum to 
Simulate the effect of the exhaust air pumps used in the 
meeeescale test facility. Primary and secondary air flows 
and exhaust plenum pressure were controlled by pneumatically 
Operated valves set remotely by differential pressure 


transmitters. 


F. INSTRUMENTATION 

A forty-eight (48) port pressure scanner, a Scanivalve, 
shown in Figure 18, (with an automatic stepping feature) 
allowed using a single pressure transducer for sensing many 
system pressures. Geofarth, {Ref. 9}, documents the logic 
and associated hardware for this system. The Scanivalve was 
employed as a computer peripheral to permit near simultaneous 
logging of system pressures. Approximate sampling of one (1) 
pressure tap/second was representative of the acquisition 
rate. The Scanivalve measured the differential pressure 
between the nominated source and a known reference. One 
Scanivalve port was open to the atmosphere and zeroed against 
an input reference signal. All other pressures were referenced 
against this port to give a precise 'gage' measurement which 


becomes a transducer output for conditioning and subsequent 
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measurement by a digital voltmeter. Pressures were sampled 
across primary and secondary orifices for mass flow calculation, 
total pressure at engine inlet, engine throat, test cell 
plenum, fifteen (15) diffuser locations and the exhaust 
plenum. Atmospheric pressure was read from an absolute 
pressure Bourdon gage and manually recorded. Pressure taps 
were sized in accordance with Reference 10. Metering orifices, 
with 8 = .7 were utilized. In order to minimize the pressure 
drop in the primary flow system, the engine nozzles were 
calibrated using the flow rate indicated by the primary flow 
orifice. After calibration the orifice was removed. 
Temperatures were measured using copper-constantan 
(Type T) thermocouples. An ice point reference was included 
in the design. Primary and secondary temperatures at 6 
diameters downstream of the orifices were recorded. Tempera- 
ture of the inlet air stream in the vicinity of the total 
pressure centerbody was also sampled. Thermocouple levels 
were input upon demand (computer controlled) to a Hewlett 
Packard 349A Seanner and relayed to a Hewlett Packard 3455 
digital voltmeter for subsequent recording. Three portable 
digital voltmeters were employed in monitoring and modifying 


the controllable parameters. 


G. DATA ACQUISITION 

An integrated automatic data acquisition system was 
employed to record fluid properties. The Hewlett Packard 
HP-IB Interface Bus under the control of a Hewlett Packard 
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9830A calculator with HP9867B Mass Storage Unit and several 
peripheral options comprised the system. A computer progran, 
Appendix H, adapted from the original work of Geopfarth 

{Ref. 9} controlled the data aquisition and storage process. 
Raw data were stored in mass memory with a hard copy backup. 
It was anticipated the data could be transferred to IBM 3033 
for processing but communication problems necessitated that 


the data be hand input into the IBM files. 
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IV. EXPERIMENTAL PROCEDURES 


Control over system operation was performed from a 
remote operating station, Figure 19. Three differential 
Peeeeime transmitters (11), (12), (13), provided positive 
control over primary air, secondary air and exhaust pressure. 
These transmitters regulated a 0-15 psig signal to three 
remotely operated valves. Dedicated pressure transducers 
provided direct reading of nozzle total pressure, cell 
pressure and exhaust pressure and were remotely monitored on 
Guyemal yoltmeters (14), (15), (16). A preliminary check 
list for system checkout and an operating guide are provided 
fieeeppendix D. Output from the scanivalve controller (17) 
could be selectively monitored as desired. Total pressure 
regulation, once the primary valve was open fully, consisted 
of remotely manipulating the compressor air bypass. 

Fach engine and diffuser combination was tested over the 
entire range of deliverable pressures as mapped against 
exhaust pressures from atmospheric to full exhauster capacity. 
A matrix of total pressure versus exhaust pressure was 
generated prior to each run to optimize the time to record 
Peeeoeand FO identify the set points for seach run. Typically, 
total nozzle pressure, PT8, was set at the prescribed value; 
exhaust plenum pressure, P14, was established; a manual code 


was input into the computer to order data acquisition. Back 
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pressure, P14, was then stepped a predetermined amount and 
the process repeated until exhauster limits were reached. 
Total pressure was then advanced and the cycle repeated. 
setting the secondary air flow to a given fraction of the 
Peemaoy air flow required an iterative process of controlling 
both flows because of their common supply. This required an 
inordinate amount of time and was not done. Instead, the 
secondary flow was incremented when desired. If additional 
data was required a dedicated run for secondary flow was 
contemplated. 

Repeatability of the data was challenged both on a 
random basis through the course of a test sequence and on 
separate dates to establish the limits of experimental 
uncertainty. Leakage checks were conducted prior to and 


during the course of each test. 
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V. DISCUSSION OF RESULTS 


As established in response to the work definition provided 
by the sponsor, the goal of the study was to design, fabricate 
and test a cold flow model of the NAPC Test Facility for 
further utilization in testing alternative diffuser geometries. 
Detailed objectives were: 

1. Model the NAPC test facility using Mach number 
Similitude and scaled geometry. 

2. Design/construct the model to allow for the greatest 
variation in test parameters. 

3. Model a representative set of engines spanning the 
Operating extremes of the actual test cells being studied. 

4. Hstablish a data base against which alternative 
geometries may be compared and provide a basis of comparison. 
5. Quantify and interpret the controlling parameters 
Which influence diffuser efficiency as a prelude to alternate 

geometry proposals. 

6. Provide a conceptual model(s) from which the second 
Phase of the study may proceed. 

eeeopecificaliy evaluate cross sectional variations in 
the straight tube diffuser to improve range and/or efficiency 
when testing the F404 engine. 

8. Explore overall systems efficiency considerations in 


the context of new design initiatives. 
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A. MODEL DESIGN/CONSTRUCTION 

The model was designed as detailed in Appendix A. The 
success of the design/construction process is measured only 
in subjective terms. The parent facility as detailed in 
Appendix C did not possess the scope of instrumentation to 
provide a characteristic mapping which would allow a direct 
comparison. The operating variables, exhaust pressure/cell 
pressure ratio and nozzle total pressure/cell pressure ratio 
Pemo@own in Figures 20 and 21, did, however, follow theory 
and closely match the general shape and bounds of model data 
Peevacd by NAPC. The full scale facility performance will 
be different from that of the model due to thermal variations, 
leakage, working gas, surface roughness and machinery support 
Seeucsure. Having satisfied Mach number and geometric 
Similitude it was reasonable to assume any substantive 
improvements in performance observed from model studies 
Should translate well to the parent facility. 

The maximum altitude achievable by the design was 
approximately 45,000 feet. The total pressure limitation of 
the Allis Chalmers was the dominant factor in this regard. 
feeewres 2la and 22 show started operation of the ejector- 
diffuser only with the TF30 and F404 in the afterburning 
mode. This altitude limitation also derives from the need 
to scale according to the largest engine. This limitation 
will obviously preclude a full determination of the useable 
feasible range of new geometries. This limitation may also 


mask some benefits of new geometries thus resulting in a 
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more conservative estimate of performance than what might 


Seeur in practice. 


B. ENGINE DESIGN 

The test engines chosen were the TF30 and F404 whose 
characteristics were noted in Appendix C. The afterburning 
mode of the TF30 was utilized as the set point for the match 
with the compressor. A top end mass flow, with the TF30 in 
A/B, of 1.863 lbm/sec was expected and a maximum of 1.75 
lbm/sec was observed. Precise measurements of the final 
nozzle diameters indicates an error of less than 4 to 1 
percent in the area ratios between planning estimates and 
the machined product. The engine design should thus provide 


over 95% coverage of the operating range of the parent 


meecrlity. 


See DATA BASE 

The 2.71 inch scaled straight ejector-diffuser was 
established as the baseline diffuser against which alterna- 
tive geometries may be contrasted. A non dimensional 
graphical representation was chosen as a preliminary method 
to interpret the test results. A gross survey of ejector- 
diffuser performance, under the influence of a parametric 
change relative to the baseline, can be readily observed. 
A detailed investigation may then be ordered to quantify any 
observable improvements in ejector-diffuser performance. 
Ideally, a real time performance map versus the baseline 


should be incorporated into the data acquisition package to 
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allow an interactive optimization during new geometry 
testing. Figures 23 and 24 are catagorized as the baseline 
for each engine tested. Improved performance will be 
evidenced by a relative displacement of any new curve verti- 
eally up and/or horizontally to the left. This equates to 
operating with higher pressure recovery for a given PT8/PS9 
which are the input specifications of any test program. The 
influence of parametric variations made during this study 
are presented in this manner for illustration. While 
conveying no additional information, an alternative repre- 
sentation of the operating characteristic by PS9/PT8 versus 


PT8/P14 is exemplified by Figure 25. 


D. PARAMETRIC VARIATIONS 
1. AefA® 

This ratio is a naturally varying parameter when 
afterburning engines are tested due to their variable exhaust 
geometry. In the F404 the ratio varies from 1.21 at IRP to 
@eeoeat maximum A/B. In the TF30, this variation ranges 
meen. 0s to 1.20. It was anticipated that, as Ae/A* 
increased for a given diffuser geometry and nozzle total to 
cell pressure ratio, (PT8/PS9), pressure recovery would 
increase. The higher Mach number at the diffuser entrance 
would govern the increase. Table 5.1 illustrates this fact 
for two runs with the F404; Figure 26 graphically conveys 


the same information. 
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Table 5.1 


Engine F404 Non A/B Bove A B 
Run No. 24 29 
PT8/PS9 dee ee e5 
Bi; PsS9 2.045 hee 


The operating envelope for any variable geometry 
engine necessitates that testing must span a broad range of 
power levels. As power is adjusted from IRP to maximum 
afterburner the exhaust to throat area ratio varies widely. 
feeemre 23, for the F404, and Figure 24, for the TF30, 
illustrate that for a fixed nozzle total to cell pressure 
ratio, the exhauster requirements decrease in response to 
better pressure recovery. The porportion, in which the 
Breescure recovery increase occurs, appears characteristic of 
Memon gine-cjector-diffuser match achieved by the design. 
The F404 full scale ejector diffuser combination shows less 
variation than the more closely matched TF30 full scale 
Sompanation. Similarly, to maintain altitude while testing 
from IRP to max A/B the exhauster must also vary it's 
operating set point to accommodate the varying demand. When 
a single test cell configuration must accommodate testing 
more than one class of engine, significant complications are 
introduced into achieving a near optimum design. Any retrofit 
of the parent facility must detail how the new geometry 


accommodates this parameter. 
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2. Secondary Flow 





Secondary flow is injected into the test cell as a 
cooling medium for the engine. The added mass saps perform- 
ance from the diffuser as a pressure recovery device. The 
diffuser entrains the additional low velocity, low energy 
flow with that of the high energy jet under complex flow 
conditions requiring greater exhauster work to sustain cell 
pressures. The postulate in the case of secondary flow is 
that, for a given nozzle total pressure and a given exhaust 
pressure, injection of secondary air increases the cell 
pressure. Secondary flow will result in a lowering of PT8/ 
PS9 or, conversely, less efficient pressure recovery. The 
experimental results are strongly supportive of this state- 
ment. As shown in Figure 27, the operating curve shifts 
lower as losses increase at the price of mass ejection. 

A detailed study of secondary effects, using the 
F404 in A/B with the 2/3 and full scale diffuser, was 
conducted as follows. Nozzle total to cell pressure ratio 
(PT8/PS9) was fixed while secondary flow was gradually 
gnceréeased. Table 5.2 for the full scale shows only minor 
Variation in pressure recovery for typical amounts of secondary 
flow. Large amounts of secondary flow have a more adverse 
Mmijeses bus GLhis is purely of academic interest as 8 percent 
secondary represents an upper bound on practical cooling 
requirements. In marked constrast, the performance of the 
F404 and the two-thirds diffuser suffers a significant 
penalty in pressure recovery. Table 5.3 and Figure 28 detail 


this observation. 
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aN Oma eRe. 


Engine F404 A/B F404 A/B 
& th, 0 8 
Run 47 46 
PT8/PS9 oS 1290 
Bi4/PS9 ao eo 
Table 5.3 

Engine Hey A/B F404 A/B 
h tt 5 0 A 
Run 30 2 
PT8/PS9 eR 275 
7 PS9 4eeas Seo! 


Expenditure of exhauster power will be required to 
achieve the same pressure in the presence of the added mass. 
A nonlinear variation in the loss of PT8/PS9 is anticipated 
due to the complex nature of mixing subsonic and supersonic 
streams. The two-thirds diffuser, having an L/D which more 
nearly matches the optimum suggested in the literature, more 
efficiently recovers pressure. This suggests that secondary 
flow effects become more prominent as the diffuser design 
becomes more efficient. The penalties in power consumption 


due to secondary flow effects are not linear, and this 
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Sleseryatilon results in wide variations in systems efficiency 
~eeotocussed in Section F. 
3. AD/A* 
In applying the Law of Continuity to the nozzle- 
cell-diffuser, a minimum diffuser area may be determined. 
The minimum area in a diffuser is specified by Ad = A* Day Pe 
where Day Pox is the stagnation pressure ratio for a shock 


x 


diffuser entrance Mach number. Allowing for an expansion to 
Mach 3.0 in the diffuser, A, (min) ranges from A* to 3.04A*. 
Matching the engine to diffuser permits upward variation in 
A, from 6.67 (in*) for the TF30 and 1.93 (in*) for the F404. 


Mmeefull scale A, is 5.768 ne which is below the minimum 


d 
for the TF30 but lower Mach numbers are experienced with this 
engine. Optimum performance for constant area diffusers, 

from Original model studies reported by NAPC, ranges from 
Ad/A* = 3.5 to 4.0. Neither of the engine extremes approaches 
this ratio with the TF30 being more closely matched while 

the F404 is undersized. As Ad/A* was varied from full scale 
femewO-tChirds, performance improved dramatically as can be 
seen in Figure 29. An Ad/A* of 6 - 7.5 appears to bound the 
gains in performance for the F404 A/B. An A, of 2.5 (in*) 

for the F404 should result in near optimal performance. No 
conclusions may be drawn for the non A/B case since improved 
performance occurs at the limit of Ad/A* tested. Static 

wall pressure profiles as shown in Figure 30 depict the 


observable changes as Ad/A* is varied from full scale to 


two-thirds for a fixed driving potential. 
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Ee F404 IMPROVEMENT 

The foregoing discussions have alluded to improvements 
in the F404 performance with variation in diffuser cross 
section, Aa An Ad/A* between 6 and 7.5 appears optimal in 
that the two-thirds and five-sixths reductions improve 
pressure recovery at all power settings. These diffusers 
can also achieve lower altitudes than the full scale, if 
that is the objective. Full scale attains 25,800 feet 
while two-thirds and five-sixths achieve 40,250 and 43,400 
feweercspectively. The two-thirds, as shown in Figure 31, 
is capable of fully started operation despite the constraints 
On driving potential observed in this test facility. The 
Gain in efficiency should be significant as previously noted 
mele ble 5.1. The exhauster can operate at higher pressures 
Mommecme same cell pressure, an obvious advantage. A ceiling 
On the potential gains cannot be ascertained from the avail- 
able data. As an example, the F404 in the non A/B mode for 
pmete/Ps9 of 6.6 would require a P14/PS9 of 1.5 for the full 
Pewee, 1.75 with the five-sixths and 2.05 for two-thirds. 
This permits a near doubling of exhaust pressure while main- 
Poroene cell pressure at test conditions. The F404 in the 
A/B mode for a near constant PT8/PS9 shows the same results. 
feeeure 30 also shows recovery occurs earlier with fewer 
losses in the two-thirds diffuser. The five-sixths and full 
peeeewattain different levels of diffusion but clearly 


greater work must be performed with the full scale diffuser. 
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In the course of the detailed investigation, both the 
two-thirds and five-sixths configurations were terminated in 
an abrupt expansion to maintain a near equivalence in L/D. 
Two additional tests were conducted with tapered afterbodies, 
Figures 9a and 9b, to capitalize on subsonic diffusion. 

Both modes of F404 operation were tested and as expected 
@eenieion is improved, as shown in Figures 32 and 33. The 
improvement at lower PT8/PS9 is barely distinguishable but 
Shows distinct gains at higher levels. Since the tests were 
conducted on different dates, precise quantification was not 
attempted. The use of some geometry to enhance subsonic 
diffusion, such as the taper afterbody, merits consideration 


Boeony retrofit proposal. 


F, SYSTEMS EFFICIENCY 

The complexity of the diffusion process makes the task 
of measuring the cost benefit of a design change a subtly 
challenging endeavor. The gains derived from a geometric 
change must be integrated over the test cycle for each engine. 
A typical jet engine test represents a non steady state 
problem where the time at a given power level becomes a 
Significant factor when evaluating power consumption costs. 
Assuming testing only at discrete power settings, the cost 
of testing at each setting can be placed on a cost/unit time 
basis and total cost summed by integrating over the time 


inmiverval ror the test. 
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The efficiency of the system includes not only the 
ejector-diffuser but must reflect the efficiency aspects of 
the exhaust heat exchanger, the exhaust control valves and 
exhausters themselves. It is postulated that only one match 
of test conditions and these system components exists. A 
shift off design as prompted by new flow conditions such as 
Memer power Or secondary flow will dramatically influence 
Geral! power consumption since it is in direct proportion 
to the individual efficiencies of each component. An 
eeiisiration, utilizing a much simplified model for the 
generalized case of testing with secondary flow and, making an 
allowance for auxiliary exhaustion of the secondary, provides 
a simple cost basing example. The test set up is as shown 
in Figure 36. An energy balance across a simple fan is 
Moeged an this case for illustration only. The total 


work done by the fan per pound of working substance is 


Hy where 
P P ee we 
2 ole 4 2 o 
H, =~ -=- —m- += +- — +Z,-24 
t oO 2 2 e 1 
0 Be Be 


Hy = 


Fan total efficiency is often expressed as the ratio of the 


work done on the gas divided by the input shaft work or: 
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f = a hC—C(<— SC constant hor Umi t 
kw consistency) 
Fan efficiency as a function of capacity follows a general 
Variation as shown in Figure 35. 

As operation shifts off design in either direction 
efficiency decreases substantially. Testing engines not 
properly matched must pay severe penalties in the cost of 
power consumption. Added mass alone provides a proportion 
increase as well. Capacity is observed to vary with the 
Speed of a fan, static pressure with speed squared and 
required power with speed cubed. 


mimewcost per lbm is equal to 


P ne 


An auxillary ejector employed solely to remove secondary 
flow must operate between cell pressure and something close 
femagmospheric. The cost per lbm for an auxiliary ejector 


would follow a similar discussion and may be described as 


ie =P 


Power (kw) atm eell 


Os : “aaj 
The combined work for the system to be more efficient must 


be less than the work of the original system without the 


auxiliary. Optimizing on a cost basis thus becomes quite 
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complex. aAs observed, with an oversized diffuser, the 

system pays little penalty in terms of pressure (PS9) for 
exhausting secondary flow. The added mass does, however, 
Pe@emawcirec, cOst from capacity considerations. A properly 
matched diffuser will cause a shift of the exhauster to an 
even less efficient setting and higher attendant costs. 


tI 
u 


in the 


oO 


Similarly, the IRP testing setting pays a lower pric 
presence of secondary flow than maximum A/B. The time 
meecor then becomes crucial to assess total cost. An effi- 
Cient auxiliary ejector could, coupled with a matched 
ejector-diffuser, markedly improve overall efficiency by 
eliminating extreme fluctuations in diffuser efficiency and 
in turn controlling the variations in the time the exhauster 
msc spend off design. 


In the absence of an auxiliary ejector, testing philos- 


jeme intervals at a test condition (i.e., IRP) are of 
sufficient duration, consideration could be given to 
meconiieuring the cell for each major power level with a 
more closely matched diffuser. This could be accomplished 
by designing a series of pre-sized liners which could be 


miserred in the full scale diffuser. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

1. The cold flow ejector diffuser model developed 
Within the context of the study provides a versatile, 
although specifically tailored test bed, upon which geometric 
Variations of the parent test facility may be experimentally 
evaluated. 

2. A complex interdependency of geometric parameters 
which influence the pressure recovery mechanism exists. New 
designs should, therefore, attempt to incorporate as many 
degrees of freedom as practicable to allow optimization of 
the pressure recovery process. 

3. When designing retrofits against a baseline model a 
real time graphical presentation of the performance curves, 
for old versus new, will enhance optimization by allowing 
the results to direct the conduct of the investigation. 

4. Substantial improvements in pressure recovery when 
testing the F404 engine can be achieved through an alteration 
of the length to diameter ratio of the constant area ejector 


Biettiser currently in use. 


B. RECOMMENDATIONS 

1. Upon successful resolution of the variable area 
diffuser vibration phenomenon, modify the test facility to 
accommodate the phenomenon and map the performance of that 


diffuser. 
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2. Using the results of the combined constant area and 
variable area studies, design, construct and test alterna- 
tive geometries. 

3. Modify the test facility by adapting the test cell 
for a separately driven ejector and evaluate in greater 
detail the added mass effect. 

4. Modify the test facility to receive its secondary 
Seeminput from an external source to preclude cross talk 
between primary and secondary flows. 

peeeexplore the possibility of including Schlerin 
photography to aid the investigative process and better 
document the geometric influences of new diffuser concepts. 
This would permit a realistic interpretation of the boundary 
layer interactions. 

6. Data acquisition must be upgraded to accommodate 
data transfer to the in-house IBM 3033. A dedicated phone 
line with modem would be the first initiative warranted. 

A real time feedback to help focus the investigation is 


strongly recommended. 
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Normal Shock 





Double Wedge With Oblique Shock 


Figure 5. Shock Strength Model 
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Pressure Taps (Variable Area Diffuser) 


Figure 14. 
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Model Test Facility Remote Operating Station 


Figure 19. 
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APPENDIX A 
DEVELOPMENT 


Design of a subscale altitude test facility to approxi- 
mate the salient features of the parent facility at the 
Naval Air Propulsion Center was governed by a multiplicity 
of interwoven factors. The underlayment for the design was 
the motive air supply; compressed air from an Allis Chaimers 
twelve-stage axial compressor (Figure 7?) eeliiversdietates of 
the air supply qualified several engines from the family of 
engines tested by NAPC as candidates for scaled testing. 

The candidate engines elected, as listed in Table B,3 were, 
from a first cut, the most likely to give a broad represen- 
tation of existing test frames suitable for comparative 
analysis with alternative ejector-diffuser geometries. Two 
afterburning engines were elected to span the operating 
range of the test facility from zero induced secondary flow 
to five (5) percent secondary flow. The choice of engines 
provided the vital ingredient upon which scaling of the 

Peer lity could proceed. 

scaling. Scaling to achieve Mach number similitude was 
elected consistent with past studies by Merkli {Ref. 5} and 
Bevilaqua and Combs {Ref. 11}. The geometry of a scale 
model may easily match the prototype but simultaneous matching 


of Mach and Reynolds numbers is impossible. A match in 
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Mach number will present a model with a smaller Reynolds 
number. A match of Reynolds number induces a higher Mach 
number in the model. Noting that large Reynolds numbers, 
Samemacent with fully turbulent flow, are characteristic of 
the prototype, any variations in Reynold number would affect 
pemenes only if a shift to less than fully turbulent flow 

was created. At a projected mass flow rate for the model of 
-5 lbm/sec, a simple calculation results in a Reynold number 
in excess of 1E6 thus relegating Reynolds effects to second 
order. It bears observation, however, that any flow phenomena 
which are sensitive to Reynolds number such as separation and 
reattachment will not result in agreement between model and 
prototype. Any improvement in diffusion which results from 

a geometric change must address this consideration. 

Once Mach number had been established as the scaling 
parameter the cold flow model carried with it a significant 
scaling bonus. Mach number will ratio out any thermal 
effects since temperature appears as a dependent variable in 
both the stream and sonic velocities which comprise the ratio. 
In the context of this study, an order of magnitude difference 
between cold flow and hot flow temperatures will fail to 
elevate Reynolds effects beyond second order. At worst, an 
error within the range of computational accuracy is antici- 
pated due to temperature extremes between model and prototype 
with the model outperforming the prototype. Work conducted 
by Welch {Ref. 16} with subsonic exhaust stack ejectors 


using Mach number scaling shows deviations of less than 1% 
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between hot and cold flow model test results. An order of 
magnitude in temperature variation occurred in these studies. 
The TF30 in the afterburning mode, having the largest 
throat area, governed the compressor-engine match. One 
dimensional isentropic nozzle flow theory for choking requires 


that mass flow obey the following expression: 


° A* x ve 
m pretty 
Vale 
0 


The available air supply had the capacity to deliver 
2.65 atmospheres and 12.0 lbm/sec at 600 degrees R. 2.65 
atmospheres would be the maximum achievable ratio of total 
pressure to exhaust pressure under atmospheric conditions 
in the nozzle exit. This ratio was below the desired test 
range but could be boosted by utilizing an exhauster to lower 
exhaust pressure at the expense of air flow to drive the 
apparatus. 

A survey of ejectors previously driven by this compressor 
revealed one design with a convergent-divergent nozzle, 
operating with half (3) an atmosphere back pressure, capable 
of pumping 2.0 lbm/sec with the exhauster drawing 8.85 lbm/sec. 
The total flow of 10.85 lbm/seec was well within the capability 
of the compressor and 2.0 lbm/sec was chosen as the design 
mass flow rate for an expected ratio of total pressure to 
exhaust pressure of 5.70. For 2.0 lbm/sec at 2.65 atmospheres 
and 600° R, a throat diameter (d*) was computed to be 1.735 


inches. Conservatively, a primary nozzle throat of 1.675 


a2 





inches was chosen, which resulted in an A* = 2.204 Gn) 
and mass flow equal to 1.863 lbm/sec. 

The TF30 has an actual throat area of 7.5 (£t7) and 
Seame ver of 3.09 ft. Dividing this by the throat of the 
me@igel, 4 Scaling factor of 22.139 was derived. Full scale 
drawings of the test cell and diffuser assemblies to be 
modeled were scaled using this factor. License was taken to 
modify supports or stiffeners to accommodate fabrication and 
assembly. Detail drawings of the scaled model are included 


as Appendix F. 
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APE RN OLA 5 
NAPC TEST FACILITY IMPROVEMENT PROGRAM 


The Naval Air Propulsion Center is a major jet engine 
imestecacility, located in Trenton, New Jersey. It is the 
only facility in the nation capable at one site of testing 
turbojet/turbofan, turpoprop/turboshaft engines under sea 
level, altitude and environmental conditions. 

Engine Testing. The engine facility is composed of 
three major divisions: the Blower Wing, Test Wing and 
Exhauster Wing. A schematic is presented as Figure 36. 

Blower Wing. they eEower Wine Ontains centrifugal air 
compressors and air conditioning systems which provide air 
to the test engine under the same conditions experienced 
byean aircraft in flight. Four 6,000 horsepower centrifugal 
blowers, one 30,000 horsepower gas turbine powered axial 
compressor, 5,000 tons of refrigeration, and an oil-fired 
indirect air heater are utilized to provide air flows up to 
700 lbm/sec, at pressures up to five atmospheres and at air 
temperatures ranging from -65°F to +650°F. With these inlet 
conditions to the engine, the center can simulate flight 
velocities up to three times the speed of sound. 

Dest Wing. The Test Wing contains eleven test cells 
and their associated control rooms. Three of these cells 
are large altitude chambers, four are small altitude chambers 


for turboprop/turboshaft/auxiliary power unit testing, two 
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are large sea level test cells, one all purpose test tunnel 
and a helicopter transmission test facility. Test cell 
capabilities are summarized in Table B.1. 

Hxhauster Wing. The Exhauster Wing contains the air 
pumping machinery required to produce low pressure in the 
altitude test cells. Fourteen of these pumps with a combined 
power of 56,000 horsepower are utilized in conjunction with 
Test Chamber exhaust ejectors to simulate altitudes up to 
100,000 feet. Table B.2 summarizes the performance parameters 
of ejector-diffuser (Figure 34) accompanies the large engine 
testing with straight tube diffusers accommodating smaller 
emeides. Iwo of the éngines which span the range of operation 
are the TF30 and the F404, whose characteristics are shown in 
Meaole B.3. 

Facility Improvement Progran. DMimgienuaeueon i982, an 
Mbeitauiye to reduce the power consumption costs, directly 
related to engine testing, was proposed. 

The stated objective was: Improve ejector-diffuser 
performance in NAPC altitude test cells to minimize exhauster 
power costs. 

The appraoch proposed was: 

Phase I. ourvey the community for current advancements 
in ejector-diffuser performance, high-temperature materials 
applications and related functional fields. Examine alter- 
nate extended variable geometry ejector-diffuser concepts 
which will provide optimum performance by accommodating engine 


Hozzle throat area variation. 
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MASS FLOW 
VabLOCTTY 


ENGINE NOZZLE EXHAUST TEMP 
ENGINE NOZZLE PRESSURE RATIO 
ENGINE NOZZLE AREA 


OPTIMUM DIFFUSER AREA TO 
ENGINE NOZZLE THROAT AREA 


500-2) 300 LB, SRC 

SUPERSONIC AT ENGINE NOZZLE 
WITH OBLIQUE SHOCKING TO 
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1000°F - 3500°F (CORE) 

a= Ag. 
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SECONDARY AIR TO PRIMARY 
AIR MASS FLOW RATIO BOGE = 25 


TEST CELL ALTITUDE PRESSURE 
SECONDARY AIR TEMPERATURE 
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Phase Il. Select one or two of the most feasible 
concepts and evaluate performance with cold flow model testing. 
Select the optimum concept and confirm mechanical and 
aerodynamic performance with hot flow model testing. Analyze 
full-scale implementation cost versus potential power savings 
and determine payback period. 

Phase IIl, Design, fabricate, install, test and evaluate 
a full-scale ejector-diffuser in one NAPC altitude test cell. 
Convert the remaining two NAPC test cells to full-scale 


ejector-diffuser. 
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APPENDIX C 
DIFFUSER PROPOSALS 


Proposals to modify the baseline diffuser geometries 
were developed with emphasis towards providing control over 
the shock mechanism. The design limitations were imposed by 
Inintaining geometric similarity of the flow paths and the 
range of engines to be tested. Whereas simplicity would be 
incorporated where feasible, no constraints were imposed on 
the design with resepcet to strength, thermal effects, vibration 
or leakage. 

Translating Wedge. A double hinged wedge in a rectangular 
duct was the first proposal considered. This assembly is 
shown in Figure 38. The two dimensional wedge was expected 
Eemprovide more positive control over the strength of the 
shock system compared to the cone centerbody. All of the 
experimenters who have investigated a second throat diffuser 
have concurred that an optimum second throat size and axial 
position relative to the nozzle exit exist. The wedge would 
allow a finer control of the size versus axial position of the 
second throat than the cone assembly. The current centerbody 
notably couples the size of the second throat with the axial 
position of the centerbody. The translating wedge provides 
uncoupling of these variables with an expectation that the 
Optimum can be approached by adjusting the second ramp to 


facilitate starting, then translating the wedge to move the 


108 





Second throat to a position of lower Mach number, which should 
improve performance. The wedge would then be mapped against 
the baseline configurations for analysis. Current design 
techniques call for running a matrix at various settings, 
shutting down, reviewing the data, developing a new matrix 
based upon judgement and repeating the cycle. Cost and time 
consumption without achieving any guarantee of an optimum are 
a natural by-product of this process. As the number of 
independent variables increases, the test matrix becomes 

much more complex with the possible permutations following 
combinatorial theory. A simplified matrix of the test 
process as shown in Table C.1 leads one to recognize the 
merit of online evaluation. A real time mapping of pressure 
ratios would be prescribed for evaluating this model. 

This would permit detailed investigations when a point of 
Significance was reached. Typically, once starting was 
confirmed, the wedge angles and/or their axial positions 
could be varied and the effect noted. 

Auxiliary Mass Ejection. The deleterious effect of 
secondary flow gives rise to the possibility of equipping 
the test cell with an auxiliary ejector. This proposal, 
while not new, has oft been dismissed as being not cost 
effective. The recent cost spiral in exhauster power 
consumption opens the topic for renewed consideration. As 
Observed in the baseline studies, the power setting of the 
engine has a dramatic effect on exhauster requirements and 


therefore, a direct bearing on power consumption costs. The 


101 





Pamie “Ce 


Variable 
Optimization Po? Pelee Size Resa) ci On 
Goal Ag 


STARTING 


Gy 3S 


STARTING 


es 


a 


RECOVERY 


A 
A 
d 


RECOVERY 


RECOVERY 


my 


RECOVERY 


ae 


RECOVERY 


= 


RECOVERY 





H = HOLD CONSTANT 
Wee= LET VARY 


A = ADJUST 


qe 





Peeeectcncy of the exhauster, when Operating at off-design 
conditions, will be less, and the blend of an efficient 
auxiliary ejector to allow the prime exhauster to function 

at or near design should enhance overall efficiency. The 
ramifications of this approach are detailed in the discussion 


Chemaecsulets . 
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APPENDIX D 
EXPERIMENTAL PROCEDURE 


System Checkout. The Allis-Chalmers compressor is 
maintained and operated by TPL personnel. Twenty minutes of 
Prelwyorication is required on the compressor prior to start 
followed by approximately twenty minutes of warmup before 
the compressor is ready to assume the load of supplying air 
to the experimental apparatus. During this time it is 
prudent to accomplish the following checks and tasks: 

1. Examine all pressure taps, tubing, and connections 
to Scanivalve port manifold and the two dedicated pressure 
transducers. Verify instrumentation is connected in 
accordance with Figure 39. 

2. Turn on thermocouple ice point reference, and examine 
all thermocouples for broken wires or loose connections. 

3. Hand test all PVC couplings for tightness and check 
to see that the primary and secondary root valves are open. 

twee lueneon the HP-9830A Calculator and printer, 
HP-9867B Mass Memory Storage Unit, Scanivalve Multiplexer 
Oo, 7 MUX), PH-3495A Scanner, HP=3455 Digital Voltmeter, 
Seanivalve control power supply, and the three separate 
digital voltmeters used for monitoring centerbody drive 
voltage, engine test cell pressure, and exhaust chamber 


pressure. 
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5. Load the program "VIBTEM" (Table i eitcoe the memory 
of the HP-9830A calculator. Run the program once to ensure 
there are no anomalous readings from any thermocouple or 
pressure tap. 

6. Read and record atmospheric pressure from the 
Wallace and Tiernan gage. 

Procedure to Conduct Data Runs. Control of the experi- 
ment is exercised at the remote operating station. (Figure 
19). *###*e2WARNING*#**##% FATLURE TO OPEN THE EXHAUST VALVE 
FIRST CAN RESULT IN OVERPRESSURIZATION OF THE SYSTEM. The 
system is brought on line by opening the exhaust valve 
fully and then the primary air may be cut into the system. 


Monitoring of total and exhaust pressure on the digital 


voltmeters allows setting of test point pressures in accordance 


with the test matrix. 
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APPENDIX E 
SECOND THROAT DIFFUSERS 


second throat ejector-diffusers have had wide acceptance 
in gas turbine engine testing due to their ability to 
provide systems flexibility to cope with the variabilities 
involved in altitude testing. A variable area second throat 
geometry such as that shown in Figure 38 was developed when 
sizing and location of the optimum second throat was loosely 
defined. The idealization of the process is well under- 
stood, as detailed by Shapiro {Ref. 8} in his discussion of 
Supersonic wind tunnels. The objective is to seek the maximun 
exhaust pressure at which the ejector-diffuser once started, 
can be maintained. A brief description of the operation 
permits an appreciation of the phenomenon involved. As 
mass flow through the nozzle is accelerated, the flow becomes 
Seeersonic and will cause a decrease in cell pressure by 
Mixing. Exhaust pressure is lowered until a minimum cell 
pressure is attained with the ejector-diffuser then being 
eomcrdered "started." At this point, the shock stands 
upstream of the secondary throat and cell pressure becomes 
independent of exhauster pressure. Exhaust pressure may then 
be increased to the point where cell pressure begins to 
rise. This establishes the system's operating range. The 


variable geometry with a conical centerbody evolved to 
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accommodate the complex mix of parameters required to 
approach even near optimum operation. This concept, while 
attractive, couples a decrease in second throat area with a 
eeieee in axial position of that Uhroat,losing a degree of 
freedom which may be exploited for further gains. Although 
the goal of the design is to alter the second throat, the 
centerbody itself will influence the character of the shock 
system and, thus, may also be in direct competition with 
second throat effects as related to pressure recovery. 
Adding a degree of freedom here may also improve performance. 
The final design of the variable diffuser utilized by 
NAPC was formulated in the early 60's, and the rationale 
behind the final geometry is not well defined. A best 
estimate is that the design was a compromise between model 
test studies and manufacturing ease and costs. The need to 
optimize the design for small percentage improvements in 


exhauster back pressure were likely secondary. 
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APPENDIX F 
SCALED DRAWINGS 


The scaled drawings in this Appendix represent the 
principal components of the design. All linear dimensions 


are in inches and angular measurements in degrees. 
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APPENDIX G 


DATA TABLES 


This appendix summarizes the reduced data collected 


during the course of this study. One set of raw data is 


included to summarize the details of the data acquisition 


process. 


The following abbreviations and units refer only 


to the data contained herein. 


EXH 
PRI 


== mead J GT) a (d=) Cd aT) se) ado] a) ef Id a 
Oo 
~ 
‘OU 


SEC 


Abbreviations and Units 
Atmospheric Pressure (in. Hg) 


Secondary Orifice Pressures Upstream (in. H20) 
Secondary Orifice Pressures Downstream (in. H20) 
Primary Orifice Pressures Upstream (in. H20) 
Primary Orifice Pressures Downstream (in. H20) 


Total Pressure - PT8 (in. H20) 
Inlet Static Pressure (in. H20) 
Cell Pressure - PS9 (in. H20) 
Nozzle Entrance Pressure (in. H20) 
Nozzle Throat Pressure (in. H20) 
Diffuser Wall Pressures (in. H20) 
Exhaust Pressure - P14 (in. H20) 
Primary Orifice Temperature (R) 


Secondary Orifice Temperature (R) 


aa 





iaecoT Total Temperature (R) 
MASS FLOW (l1bm/sec) 
Pao lAG Gi neteabs.) 


DeoTAG (Degrees R) 
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APPENDIX H 
DATA ACQUISITION PROGRAM 


A computer program which details the data acquisition 
Beeeess is included in this Appendix. VIBTEM was executed 


on a Hewlett Packard 9830 and is written in BASIC. 


te 





NASM Lag 
“SHO 40 


a 
bre bk 


sda. 


HI 





CaP Je (CSP Is *C pede Ld. 

PuHltd dO YH34* SHO HLHOW 45LH3., 4310 

aa2=2 LHW 

A32=H LYb 

432s LHW 

43é=), LYW 

da2=% LHW 

COS IS*C BS IWC Stl ae) sas 1 WI 

LHIdd 

LHIdd 

LHIdd 

LHId4d 

LHIdd 

LHId4d 

T3aHHHHO Hda=a WA 

SJHIQHS4 WAd= W34 

VANHYHD WALLINSHHaAL=3 W3e 

TSHHUHO THe2o Wad 

VaHHUHI O1=13 Wad 

# ASHWHIS=5 Wa 

UI adel a dAHHBIS=$5 Wd 

"HOTLI9S AanLY43dW3l 404 SASYTaH Wad 
aes dats = $6 W34 

140d &/S LHASSdd = d Wad 

LdQd HOTH = 2H Wad 

Laod Wo = TH Wad 

47S W341saq = 4 Wad 
“HOTLISS 8/5 404d SIGHT aus W3d 
Wad 
eT-T STAHHYHS W3d 
Ja HO SLHAMSAN SHAM Aa Lea aa SWaO sda OSH LI Wad 
oP= | ee JAH Lad HASMLAd 6S_ SA THAIHHIS 40 Wad 
SUIHHHIS WE LHSN3S SHAOSaS4 WHAM STH M34 
*HOTLd1 49530 Wad 
Wad 


wHIALITA, 83 


WYH 311s 


Wad 


on 


=) 
Zi ™ 


iene uo a 


oO 


a 


Natt Fa 





SoS 
OC ee 


ce ra oy 


Don Toe 


Va Ca Op coe 


ay ete a 


aj cn 7T UD 4a 0 0 ee 


oy 


5 owe 


ned 


—, 
RD 


To iT 


tit 


Wa tt et et et et ee 


aed 
at 


eT 
mal? 


ate ae Fe 


DD WMO So oo 


Sogo uo ree 


oT Ee Tee re te Te 


134 





+s, 


See ER + ee eee Hie esse eres 


IRE 


i] 
PRE 


os 


mule 


PET 


F400 


F 





~ 


ba 





By 2x 
HE: 
OPT Hoos P a. 
SEGMENT To 


Ly 
IF H 


we 
° 
a Ci 


le 


REM#*4422START 


WRITE 


= 
Los 
os i) 
eed rt pris 
Se cae or. me i 
w. II = a} hea} er 
4 fe oo 
Ute i 4 Ta Dt wr = 
or a Ese <f oTi Lil ae 
is = — Ls — Ess 
Sl _ = Bcatad i 13 
mn se asl [BE nm 8 —_ = 
seat on We Lil einen (ee 
+ = = Lil ck P-L (ms th 
wT te B= om Tos 
ww oe tes = cen St ie 
Bocuse ae ties Lud = Li wat me 
1 aa Lu bem . kK T+ 
iT ce te lil cd Ts i a 
co Lut bij = T _ att 1) mee 
fo eae camera Dem] Morte mn Lad 
a lu f m0 to ° = = w= om lu == 
a) Et Reet ees man ke — UD 9 Ce oe e 
wae FS Tw coc =P Le wt OC ee 
~ Lu Bo eo 3: i = eet 2 

-— = Ih = bem be be eR et — i 
Lj lb b— em rere E CC cr = | m— Hi LU oo Re he et 
capper) UP (6 eet ono aes eee LD eb ed Rm fame eters) ce aoe 
Sees Tot hee EE eer re ee Re ee — 
LEDs ad SS Take Pad eames cee eet eee TOR LES Leb merit Rentl lize!) Mee Vet ey UL est geg adits Semele fee Soe LS Lee eS 
SuwWki aA Ss ae ahaa ea KH CSL eek 
ete Ta ee Te Ty ve eee ye ee Te Te TY re ee ee ee Te ace a ZT 
Pom eT TT ee ep ge bate Pe ree ee Te ae ea oh wee eg eg eg eg ce eh 
CC 0 <P oer poop osr oe owr och oer oer uu 0 te a ae 


1) 


3. 


tele 





t 
oj 





" PORT’ s 


so 


Wye 


#" SF 


AHITVALE 


1 7 





r | 
= 09 UD oD 
-. etwts 

a 

ee = 

my Setter) tikes 

T ae rea 

Hey et ee he 

PAP echo ae tt 
; — = 
eae 








A 


=H 


Hi 


IF 





MAT 


Fk 


Es 


SE 








Be 


ILTAG 


¢ iH oe) 


PEHTER TIE 


i] 
LHFLIT 


910 


136 


moe 
t 


‘HES 


sles "TH 


ee 


IH 


‘TTIOC 


Tee ee oe 


1 


HOF VTE 


“EHTER! 


ee 
Pe [a | 
— 1 
ges 
oe 
eee lient) 


REFERIT “s 


1) 


ips 
iT. 


F 


oss 


Be 


ME 


] 


q 


4 


4 


pee) EA se 


Tt Ts w— 








+> 
ar i aa 
Ze ia ain) oe 
4 a 
_ uD 
my 
lid til = 
ale a mens 
Ree IM 
we 
Mm ml 
we os oe uJ 
Ao - 
roe 
La Le tu Lid 


iy SD OM aE es 
f-— 10 oT. 
yy i iy wt 
ot wet ot et 


ara 
iJ 
es 
= 
a 0) me 
“. [a 
wo 
lu t =e 
af he =s- = — 
LJ Lames. © wn —. 
= mee y= un 
7 ry My — st —s"~ 
4 Ses um sk rc woot fe 
n maa nw an iu COM nm 
Ca Wwe oP mus wea ak 
im eT 4 wt 1 Tr ee bu ow IE ee 
wl wt es m (Ske Ae + 
TS Wi [1] ome we WD ane 
tal wm Lud tL > rer + 
Ee eee eE OD te 1 ee 
Hee Se eee EEG Owe o- 
SCOWOee CwWowWwSE ot 
ee eae Ise heen aca eee eae eed ES ete od Pome eee Lee 


aT Me eT ae oT ae TD OSD Te Tae TD Te ED aT Ee ws 


MPC Wap wp Pe wm eye ap ote. 0 


et et et wet we wee EE ee ey ey ey 


Se) eS) Se 9 oS a a Ss SS aaa 


Wey 


EEE SSS ERE ERE EEE ER Ree eee REE Ee 


ow 

lis 

— 

ra 

= 

3 

iu 

a 

an ™ 

ua - 
mr 

~ lu 

Lo ae 

u4 a 

Wed Fé 

_™ saat 

ic 

ie wad 

ua co 

n ed 

c? cre 

= Ro 

Tae te 

om 

Lig te 

Fat 

me 

CK iu iW 

ii 


yoo oF 





’ 


‘ 


, 


my 
= 
un 
tile 
ig 
f= 
0 
= 
cy 
fa 
weer 
rey 
ot f-. 
UD oe 
Col oe 
au 
Die) ws 
. Cl om 
m™ 4 oo oo 4s 
Cou Oe Se S 
ms I “= 
+e ke be bem tT a 
wr rer ell ame wee ae 
Sle ee iif. s 
Lee rereoarun Wi 
= IO Dus Hm] He 
AunWwh EO Moo 
SE) EZ) le Ty as 
reg ee mp Te ae Pm ale oe i oe 
re ch CA OO op mr of 
= 4? wt 





+ if} 


uJ 
’ 









+: 
S Ts 
te m 
Cal} 
. uo 
4 uw 
Lil re 
i] me 
Th 
ct Tat 
l= == 
s or 
e. = 
a] Sod 
ow —_ 
2 r 
= re 
ix _ 
oak ue: 
hans Sol 
: rag 
“ wy 
en (Zak 
oy a 
n~ Da] 
ic ae 
-. mu Ce 
~ way ey 
- =: a= Ii 
= on. * <a 8 
. wed are) it 
re] 4 { am if 
iL un = iD umn 
o w. ie a ir 
S ns Tr Ite) es Wt 
tt nm + — =< wt HT 
if way uy LL | w— 
ce uu a =i Ze Linh 
tu 4 7] ian To] wo yh oe LT 
Ces os “. i 4 Ze me wT ime. 
pate age my Dice } mt coy =" ow. ™ ome go 
er oy oo +r wee gle Ute ey Re Re z= =e Te re 
Thats cA rs Fr Ty nw jem =P WT cc ify 
mT ut uo —_ om —_ oO mm oT eigei ir tT fis 
ws wr. nm Ba meds LY wt Oe —_ + it. + 
ao af co om ”~ Mee rerum ew cm, a om [ro 08) fr 
usc ee ee Ree EO TIMt We Tee eH 
=< lt = we SE ee we tet me pa me we 1] w=. Ss we fe LD fe 
‘ 3 re Ljw mt wt OT os ted er ” be iD bie ih 
me it wee or - Wo Te Lowe HH Doe wt Be bh e+ ke He oes eer 
Ww mes ls ce ust ft os er ee TP oo Ww oe Le oe ce 
= a a EE ee ees VY YT pee TV ee Ie Sue ro wa 08 Do a 
SE ERR Ae Ae OS eae _ mem Oost ot eo 
rmIOTMI SEBS Tl oWeBowy bw kK DW ke ODE We Wee 
SuUWOUODUOD UW Dea err eRe RK Sa TSE VE Awe 
WE) Ta ap TE ep a TE 
Moet gf. Ch ee eG ep oe Oe eee ae ee ee Oo a 
spf er oer toh 0 6 ou ee ae ae a a a a Re Pe 
ee et ee we ee ee te et 


3 





nee 


Riles 


ee 


LIST OF REFERENCES 


jew te or., Lustwerk, H., “Supersonic Diffusers 


For Wind Tunnels," Journal of Applied Mechanics, Vol. 
Hicmmiee 2, pp. 195-202, June 1949. 


Lukasiewicz, J. "Diffusers for Supersonic Wind Tunnels," 
Journal of the Aeronautical Sciences, September 1953, 
Boemo ly -o260. 


Hastings, S.M. and Roberts, R.C., Analysis of the 
Performance of a Two Dimensional, Variable Area Super- 
sonic Wind Tunnel Diffuser with and without Scavenge 
scoop and Model, NAVAIR Report 4384, White Oak Maryland, 
May 1957. 


Panesci, J.H., German, R.C., Improved Methods for 
Determining Second Throat Diffuser Performance of 
Zero Secondary Flow Biector Systems, AMDC-TR-65-124, 
July 1965. 


Merkli, P.#., "Pressure Recovery in Rectangular Constant 
Piesaeoupersonic Diffusers," AITAA Journal Vol. 14, No. 2, 
September 1974, pp. 168-172. 


(rms ei sds, Bitlie, Hoo. "otructure of Shock Waves! 
Pimcvlingrical Ducta," ALAA Jeurnal VYol. II, No. 1- 
Meueber 1973, pp. 1404-1403. 


Cimon, Jal.e, Supersonic miectors, AGARDograph No. 163, 
November 1972. 


Shapiro, A.H., The Dynamics and Thermodynamics of 
Compressible Fluid Flow, Ronald Press, New York 1953. 


Geopfarth, R.N., Development of a Device for the 
Incorporation of Multiple Scanivalves into a Computer- 
Controlled Data System, M.S., Thesis Naval Postgraduate 
peneoly Monterey, CA, March 1979. 


American Society of Mechanical Engineers Supplement on 
Instruments and Apparatus, Flow Measurement, 1959. 


Bevmlgquayerat., Gombs, Cora, [heory and Practice of 
Ejector Scaling, Paper at Ejector Workshop for Aerospace 


Applications Dayton, Ohio, August 1981, (ADP000530) 
Martian, B.W., Baker, P.J., "Experiments on a Supersonic 


Parallel Diffuser," Journal Mechanical Engineerin 
Seiences, March 1953, pp. 169-174. 


re 





McLafferty, G., "Theoretical Pressure Recovery Through a 
Nomnaeonoce in a Duct with Initial Boundary Layer," 
Journal of the Aeronautical Sciences, March 1953, pp. 
169-174. 


Moauaects, Voie, Dutton, J.C., Addy, A.L., "Pressure 
Recovery in a Constant Area Two Stream Supersonic 
Heeeicer, ALAA Journal Vol. 20, No. 9, September 1982, 
pe 1308-1309. 


Johnson, J. A., Wu, B.J.C., Pressure Recovery and Related 
Properties in Supersonic Diffusers, AFOSR ipo on oO 
COo2 Apri 1974. 


Welch, D.R., Hot Flow Testing of Multiple Nozzle Exhaust 


PoueuOreoyouliio, Meu. Foects, Naval Postgraduate School, 
Monterey, September 1978. 


140 





CVE eA tepiOTRIBUTION LIST 


Defense Technical Information Center 


Cameron Station 
Alexandria, VA 22314 


Library, Code 0142 
Naval Postgraduate School 
Monterey, CA 93943 


Department Chairman, Code 69 
Department of Mechanical Engineering 


Naval Postgraduate School 
Monterey, CA 93943 


Professor Paul F. Pucci (Code 69Pc) 
Department of Mechanical Engineering 


Naval Postgraduate School 
Monterey, CA 93943 


Assoc. Professor R.P. Schreeve (Code67Sf) 


Department of Aeronautics 
Naval Postgraduate School 
Monterey, CA 93943 


Mr. J. E. Hammer, Code 67 
Department of Aeronautics 
Naval Postgraduate School 
Monterey, CA 93943 


Commanding Officer 


Aiwemunons Mr. Vito Trugiica 


Naval Air Propulsion Center 
Trenton, New Jersey 08628 


CDR. James W. Molloy, USN 


1253 Elm Road 
Eueramore, MD 21227 


GDn thomas H. Waish, USCGG 


1908 Carriage Drive 
Walnut Creek, CA 94598 


141 


No. Copies 




















Ue oe 

Let RL ore 
Yer Wea patete + 
bara Ler) 
Rae Soha he! : 
4 And ostet 


eceeeeemees § Scaled ejector-diffusers for jet engine 
oy 


by {i HIRT vee q HAHAH mea mm Ha: A } 
| HII Ht | | ill l| | | 
Ne HHH) WTI MK | ; we 
ATT HIN HI Nt! i Hl I A an 










































Per est rts 
mT STW Rare 
Srey ie sy ae 

4 ili bea ah 
Tedeae bole on 
et ne ware er 

PETES T pore ty 
ik pied gE ame Mae B.S 3 wh es 
BAT OTe aw ren oer a ra 
eS Ley) OLAS Res Se RR ra 
inte Pritt he pee re he iT em Tee era aa ee eaerd 
Eee hee TS Fetter R RC RE PEI wT) 

ig er Ral UP SE ay Re oe any cree} Ce) 

Ani Set TLL eee ey Ee Y CRD nf 
pt gn eked chat te Lnes Pe ae ery pr rer 


PST SEES OREN I ERY RE aD aa 3 2768 002 04677 3 ta 
sary SE aN SEPP DY apy: DUDLEY KNOX LIBRARY 










ag wlhe ene. 
uber walrt heehee Lae 
parla deter ek dotal 
Poa pee sere 
ae Tot 










Ey 
Se me ee 






































Sel da ee abled due adele Lies 
Rae ademas wees eee 
Meira hele aaa tL Pe) 
pebieabpul iat kadai keke Tuk RY 


aa 











ong Pr eer 
dt Bulb ap ed at yt 










a Pe 











errr re 5 5 O 
BRE eee See bic Yow eet ahi yee recta f ss A 

Pare CY are ey ene Pea ann raters Pe pet A / 
Le apes etek tok LT Cit eee re 





Rear er soe 


Tae TNL CL ee ra eee roe 


a 
SY bebe danse tania aA Ie me et rd Bob tee 2 
OE EPIC RO rem tee ir seah Ta gan 
Mh pi EA BD OVROM Dm RAE beg reme = CPO ie canes 
3 Surpeeiow ah oet dub ene TTR TT ee Tere phi eee 
Pe Tat Deion TY nr ea os I) ih cet daybed 
ot pw idepereeletuadlde hare beniehyd wt Sera RE a ee i 
ula etl a ae SRG ICL PT Oe Pre 
Spreader te Ser srl rt ary 
onehdinde tee Ler ee ern 
eb nee eri ER LE ERT Ee 
Pima pepe leet PERT N S/O Ree mare aero ade 
epee epitngin ahd] tote ae Oe int or rea aaa ary ry 
List nipaet gh ata ee es 
ciate is 
aR hee meee IDI my Eire arr tabat 
Veale wpe oe EE re ee ‘a 
Lee eteder ed icine ee ee ver Ty 
oe Scare thei k fe ey RE LET 
cumini del ated LL ere i 
eee rtd ees 
Feb plke tacec pated foe et nT Pay peel 
Poel ta epel Pelee eh meer pate eo irr el 






















































Pern 

























Sp abil rad beetle ee Eee 

A dntw Dowd lite ee 
whectyion tate Series eater eT TRIE) 

Rye oade meinen eaters as bases 
ene 





CS ae ney Were en 
ot amen erie Emory ws 
Sanita EL rr Pett 
SEY at eee ea rae 
eae tare aie Eien] 
Sap biel ter sated 

















Ceres oe reir eae a 
in ead ire eee 
: 


ee 
ee 
Pere 



















Parse rie 
ene tes 
Pr Pere aT ors say 
PRAIA ih wy ny Re 
a ein ot ee 
LIS eee: 
=: Papo ye fk, 
re 













ei 


ete 






















Se Ug 
Pelt A ne ¢ " 
bg arian : - 


vier ts . 


tials 

eat Petter 
Caron 

eA t 





ttiedehite th acy Ch | 
Pad Hiatt hls ira igh att, te 
oo teat oe eae 





gab taal he Why lin anak Geni a 

ee etre ris at te Hepa 

SON Pte Ake 

mont Nady asa ey ee 

~ oat picla htaceten Reba 

Oh eer ty le eee aah 

hy Salalah atelat Glan te tah It 

negate hie nfemtyik a btca a’ 

Se A 

erase eT 
os, apag ty Yee intake 

Seip tat Ppr We uKb Tee era) 























ke Bathe 













‘ Sy lel haa Ce at ee 
Jaye thet Rade aA OOS Gar eet Ned ee ey aA Wet Seb d 
RT tit ep mak Eiiptt eade AS 
eee nae veer hy ad etek Diya wits 
ve 








“ ata r et 








Cr Seis by 
Se A Sok acai ler a td Perret 
yea ABORT eS 






re 









we Li 
Recent 
4 a Sha Vale! vi 
Pry cease? 
rege NN NE 

Siebel de "wy 
Roose NLDA 
Sena rre nN 
UO 
































Spats 
ey x 
SX tae th 
Reh 





+ 
Sneak wade 





Ot Tee PE 





KP) 











